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ABSTRACT 


A Droblem in the linear theory of elasticity is con- 
sidered wherein a layer with a circular cylindrical hole is 
subjected to a nonuniform radial displacement. The defor- 
mation is imnosed on the cyl indrical boundary such that axi- 
symmetric disolacements and stresses result. . The solution 
utilizes Navier's equations of elasticity. These equations 
are solved by use of extended Hankel transforms to obtain 
displacements. Shear and longitudinal stresses are obtained 
by transformed stress-strain relationships. Radial and 
circumferential stresses, however, are obtained directly by 
use of stress-strain equations. 

The solution o^ a problem where the imposed radial dis- 
placement is a linear function of the axial coordinate is pre- 
sented. Numerical results are given in granhical form for 
two different ratios of hole radius to layer thickness. The 
infinite integrals of the inversion formulas were evaluated 
numerically using Longman's technique for computing infinite 
integrals of oscillatory functions. 
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CHAPTER X 


INTRODUCTION 


In 1964, Scott and Miklowit z [1] * introduced the use 
of extended Hank el transforms in a problem of elastic wave 
pro po gat ion. The problem was an infinite plate with a cir- 
cular cylindrical hole. The cylindrical boundary was sub- 
jected to a uniform step radial displacement and the result- 
ing axially-symmetric compressional waves were studied. No 
attempt was made to determine the state of stress in the body. 
Seco[2] introduced the use of this type of transformation in 
solving an axisymmetric heat conduction problem in 1969* 

The extended Hankel transformations employed in this 
thesis are based on an expansion formula discovered by 
Weber [3] in I 873 . Orr [4] rediscovered Weber's formula in 
1909 by a method of contour integration. The formal proof 
of the Weber-Orr expansion formula was established in 1922 
by Titchmarsh[5] who, in his book [6] , broadened its use. The 
following expansion formulas are given by Titchmarsh[6] s 


* I 


s £• fans.} «s/s 




•rftfr) £* Skn*) <s s/jr* 









( 1 . i ) 


( 1 . 2 ) 


^Numbers in brackets refer to the Bibliography at 
the end of the thesis. 
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where 

JB& S J* tfff"') y, - J, (sm.) y» fsr°) 

& fer) 4.) « j, Cs-r) ¥* €s®) y, isr) . 


The same set of transforms used by Scott and Miklowitz 
[l] are employed in this thesis* Defining the zero and first 
order transforms as 


(£.[/<'/-)] - /if*) * £ r/Wi *• fr,V*i> •ff (1.3) 

<S, [//r-S] - - C* (!•“) 


respectively, it follows from Eqs. (1.1) and (1.2) that the 
inverse transforms are 


«•' m ■ n &fzs* ■ 


(1.5) 


( 1 . 6 ) 


Properties of the transforms which are applicable to 
this thesis are presented in Appendix A. 

The problem to be studied in this thesis is an out- 
growth of the work done by Scott and Miklowitz. The aim of 
this study is to ore sent a method for determining the static 
displacements and stresses throughout the elastic body posed 
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by Scott and Miklowitz for the case of a nonuniform radial 
displacement® 

The statement of the problem and its analytical solu- 
tion is presented in Chapter II, The technique used to eval- 
uate the inverse transforms numerically is presented in 
Chapter III. Numerical results are given in Chapter IV and 
are followed by the Summary and Conclusions in Chapter V, 



CHAPTER II 


THE PROBLEM 

Consider an elastic layer of thickness ^ with a cir- 
cular cylindrical hole of radius <£. The cylindrical coor- 
dinate system 7^0 jZ is used, T* being the radial coordinate 
and 2. being the coordinate along the axis of symmetry, as 
shown in Fig. 1. The layer is resting on a rigid foundation 
in such a way that zero shear stress and zero normal dis- 
placement exists at the interface. 

An axisymmetric radial displacement is imposed on 
the cylindrical boundary which varies . linearly in the in- 
direction. This imposed deformation, shown in Fig. 2, gives 
rise to axially-symmetric displacements and stresses through- 
out the layer. 

The problem shown in Fig. 3 represents an infinite 
free plate of thickness zA with a circular cylindrical hole. 

The plate has an axisymmetric radial displacement imposed on 
the cylindrical boundary which is also symmetric about its 
mid-plane. This type of problem, which could be of practical 
importance, is completely equivalent to the problem represented 
in Figs. 1 and 2* 

The Navier Equations. 

The Navier equations of elasticity in cylindrical 
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figure 1® The elastic layer 

geometry and coordinate system® 




Figure 2, The imposed radial 
deformation of the layer® 




Figure 3» 


The equivalent problem 



coordinates for the case of axial symmetry are (see, for 
example, Reference ?)* 


8 


- -fi \ r |~5' + - O (2.1) 

- *M jh (2.2) 


where a and are radial and longitudinal displacements, 
respectively, and y is Poisson's ratio. 

Solution by Extended Hankel Transforms, 

Following Scott and Miklowitz 1 s work, Eqs. (2.1) and 
(2.2) will be solved by employing the (jE^and transforms 
defined by Eqs. (1.3) and (1.5) and Eqs. (1.4) and (1.6), 
respectively. Using the properties of the transforms given 
in Appendix A and applying the (2, transform to Eq. (2.1) and 
the (^transform to Eq. (2.2), the following differential 
equations results 

-•<2 ( /— jt /) s*" & 7 s - & <7 — S’ QSS *ss (2.3) 




S*U> 


/—<&*) £^60 f - 3 * ££ 








(2.4) 


where "D S, 



From Fig. 2, the imposed radial displacement is eiven 
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by the relation 
Therefore, 

^ u (&,%.'$ 


(r-b) 


3£ 


s®, — TM’M 


( 2 , 6 ) 


The requirement that the cylindrical boundary be free 
of shear stress implies that 


^ w 






(2.7) 


Substituting Eqs. ( 2. 5) *02.6) and (2,?) into Eqs. (2.3) 
and (2.4-), one obtains % 

— g s A ' £F ~ sj>ia> 




( 2 . 8 ) 


f's-zp) s*m? 


57>W ^ 


(2.9) 


It can be shovm that these equations are equivalent 
to the following fourth order, uncoupled differential equa- 
tions * 

(t>*- «iVf) u . spL i * 


2 s V •Kr A w • 


/r- 


773W 


(2*10) 


( 2.11 ) 
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The general solutions of Eqs. (2.10) anil (2.11) are: 


&& Z e +■ 


Cy ^ •/" 


t- 


Substituting Eqs. (2.12) and (2.13) into Eq. (2.8) 
gives the following relationships between the arbitrary 
constants s 


f'g- <cy 




Equation (2.12) now becomes 


£T S - 5C, 






(2.13) 



Equations for Stresses* 


The stress-strain relations for axial symmetry are 


for example. Reference 8) 

S 




(2.19) 

<T^ - (X+$u) ~z 

- tt) 

(2.20) 

Oi - (X+&U.) ~ 


(2.21 ) 

<** - ✓“/ -5T 

* T%] , 

(2.22) 


where X and /U. are the Lame ' constants. 

Applying the ^©transform to Eqs. (2.19). (2.20) and 
(2.21) and the transform to Eq. (2.22), one obtains 

T>vS3 \ (•£$ u (a,*) +. s IX) ( 2 . 23 ) 


( A jT '—J. *> y-3tT- | ' + 

+ Xj^Dio 4 - torn*) £$ Jr-J ( 2 . 24 ) 

(X [ f *&s*) *• X * 4 *) - 

■ ^ A */rJ 


(2,25) 
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Equations (2.24) and (2.25) are not in a useful form 
and, therefore, Eqs. (2.20) and (2.21 ) do not possess suitable 
transforms. Radial and circumferential stresses must be 
obtained from a method other than direct transformation of 
their respective stress-strain relations. Using the integral 
expressions for Lt/r^-z^ and and Eqs. ( 2 . 20 ) and 

( 2 . 21 ), one obtains! 


Q7- 





+■ 


\ I / S_ 

+ y ” j,%f) +-y,*&s) 

. X 4—( — & ^ (2.2?) 

4 */&*) »*■ )6 Y»t\ 


S\ \ / ( * ** Sx,z') 


Ws) 




£ 


-h ySfrs) 


\ * 1 

/ s tu /r, as.) a/s 


^ J* /as) y & /as) 


( 2 . 2 B) 


If the following three conditions are satisfied, the 


order of integration and differentiation may be interchanged 
(see, for example, Reference 9)1 


(i) the integrands are continuous 

(ii) the given integral exists 

(iii) the resulting integral is uniformly convergent. 
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It can be shown that these conditions are satisfied. 
The first two conditions are obvious and the third condition 
requires knowledge of the arbitrary constants in Eqs. (2.13) 
and (2.18) which are determined in a later section. 

Interchanging the order of integration and differenti- 

/ 

ation and performing the necessary operations, Eqs. ( 2 . 27 ) 
and (2.28) become i 


^ f v /x . ( m A cSs 

+ J,*f ar )+y t 




7*j) i-ysY**) 


2^ 


J'fa) *~y, x ( 


e*S__ 


rr- _ \ (~£ + SI fs.z.Vl & 

01 ~ 


<Ss 


4 - 


+ £J±Ifd7l El 

" r * J» +yA 


y*'&<> ) 


where = s' [ *T & + 

A. (sS) = e‘*[ -sz) C 3 - 


(2 


(2 


(2 


S £*, J -f € - sz) c* 



. 79 ) 

.30) 
• 3D 


(2.32) 
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Completing the operations as prescribed by Eqs. (2.23) 
and (2.26), the following expressions are obtained for trans- 
formed longitudinal and shear stresses? 

- X - S' e SJ£ [ f 2- &2J 


S ( C,J 


' SZ [ f 2 - £>V 




Y SZ (l-ZV)^ C4. *• S C l~2.it) C 2*] ^ ( 2 . i 3 ) 

<5^ = 2yU ^ e 3Z [ sc, ~ - $z)c 3 ] - 

- e.- S3L £ sc z v- (/-£*> y-s-z) cr^J 1 2. v‘+) 


The longitudinal and shear stresses can now be obtained 
by inverting Eqs. (2.33) and (2. 3^)* Using Eqs. (1.5) and 
(1.6) » one has 1 


(rrr) 



S ( S e !L) tat m a/s 


J, ■t-y/fax) 


^~r% 


(rzz) - \ 


S S^aCXz) R. t els 

J -t- . 


(2.35) 

(2.36) 
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Determination of Arbitrary Constants . 

From the statement of the problem, the following 
boundary conditions can be deduced (refer to Fig* 1 ) » 


ci (7703 = 0 

(2.3?) 

(t; 0*3 =• 0 

(2.38) 

<&x (n *3 - 0 

(2.39) 

cv (-r, A) « 0 . 

(2.40) 


Applying the (^transform to both sides of Eqs. (2*37) 
and (2.40) and the < 2 » transform tc both sides of Eqs. (2.38) 
and (2.39), one obtains for the transformed boundary conditions! 

( 5 , 0 } “ O (2.41) 



(S.o'i = 0 

(2.42) 


- O 

(2,43) 


Ur<Ts,d) * o , (2.44) 

Substituting these boundary conditions into Eqs. (2.18), 
(2.33) and (2.34) , the following equations in C ljf Cg C 3 and Cq, 
results 
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-sc, - sc* + zo-&) £3 - z &-») c h 


= O (2.45) 


S'C, - SC & - fj-£&) C 3 - (/-gif') Cy, - O 

se s ^c, - se^ - ^-^sJi)e s ^c 3 -{/-&+s*L)e*^ 
« ^ ^ ~ ' — 5 : e ^3 SF e c f s 


(2.46) 


*0 ( 2 . 47 ) 


/rs 


s(2.4d) 


Equations (2.45) through (2.48) have the solutions 


c, “ ■ 


(2. 49) 

rrs 9 C - «*** +■ *-***) 

<r* - 

739 *$ f f /-**>) -sAf5~yy)J A- e^fSsnty) HsriJ J 

( 2.50) 

/rs* s’/-y)('- Y*A - 

^3 - 

7**/ <A f ~ J 

ats* <r*-y)f-4s4-e *•* ^ 

(2*81) 


7?*o & / e ** - r 

(2.82) 

V 

srs* ?. e^A) 

Special 

Case , ^ *0 . 



To reduce this problem to that worked by Scott and 
Miklowitz, i.e.» one where the imposed radial displacement 
is constant with respect to Z. is merely set equal to zero. 
This causes the arbitrary constants to vanish and Eqs. (2.13) » 
(2.18), (2,29), (2.30), (2.33) and (2.34) becomes 

z A 


u f yjs } - 


JTS 


(2.53) 
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^?VjaO * 

O 

5 

® 

CM 

Q^Cr^z) * 

i J, 8 ^) + 

(2.55) 

ft 2^ « 

3LM (** s ZCs^ & Jr 

^ & 3 * ^ */2w^ 

(2.56) 

of 

* o 

(2.57) 

°rz (rjs) 

m. O , 

(2.58) 


Using the property from Appendix A that 

£ 


& __ 
L S £ “ 


— *** 


the above equations can be summarized as« 




(2.59) 

£*/*“) * 

i*g, 

7^ 

(2,60) 

o^» ^V~) ~ 

r* 

( 2 , 61 ) 




6 ^ 


( 2 . 62 ) 
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The expressions for U and W> agree with the plotted 
data ("long-time solution") by Scott and Miklowit z. In 
addition, the expressions for displacements and stresses 
agree exactly with Lame's solution for a thick-walled cy- 
linder (see, for example. Reference 10) when the outer radius 
is infinite. 

Nondimensionalized Equations. 

The equations for stresses and displacements can be 
nondimensionalized by introducing the following dimensionless 
variables s tz( = «A, //- , y . sA, />-- 

With these substitutions, Eqs. ( 2. 5) » (2.6), ( 2. 35) » (2.36), 
(2.29) and ( 2 . 30 ) take the nondimensional forms t 



Ufa f) 

L JSfar) fr?) 

(2.63 ; 

. 


(2.64) 

OS Cfif) 

£ 

Y 5f C*LAV) 

(2.65) 


r“ 7 

(2.66) 
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^ J* J'Y**) ^y/Yvf) 

■*£ ( ** £ &f) £» C*t 4 /l#) e/^ 

/ ^^7 

/O c*~**) 

TO ^V;») _ *> f V^Vyo XiQS.tj/,*)'* 

£■ ~ <-'*»Xhp)). j,'fy r )+ y, '<wr) 

+ iJ^Z , 

where 

- C, e fr *-c t e e *r c s fe rr +■ 


( 2 . 67 ) 


(2.68) 


4^?0 . f e&lz. . c] e* + 

+[ CZ&3-+ *]e* - 


5 ^* 





8 

^yS 


CSz ^ ^[ +**'*- fMt*t)f?)C s - 


&*J - e *£ fi £ 




#* 



i 
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^ (f,y) * # fr £yc, - o-**-?*) 07 - 

- e ^ r f ycg c ¥ J 

c « '^AJL + £l£l£^2si2. K 

/ ft-tt)f- #*?- e ^ e"**) 

£* _ f e *£- O-**) /- tf^CSHns) rfj £ 

* ~~~ f~ g.-rry— 

c — dz^L e ~*% 

5 /r? & * 7 - e**+c ~*f ) 

CT" s ^W// C+S**?) “ g^jt 

* try* f/-#y- 7 l 7 -e*? * e-**) 

fi&V) * y'fe^fa rfr*) +■ + 

A {&¥) * 7fe fr [f*0-&S)-fr)c,--yc,j + 

+ &~ f? [f'-zS'-**>) -f<?)<f -yc t JJ 


^ ^ * J* <??) K <Vr) - J» <v?) y#4^) 

F (<av) * \ {/>?)% fag) ~ ^*(«ty)£ff*i) . 
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It can be shown that a singular point occurs at 
* A » 2 * • Radial and longitudinal displacements 

and shear stress have imposed boundary conditions at this 
point* however, radial, longitudinal and circumferential 
stresses increase without bound. 

For large arguments of the independent variable S * 
the functions p and pg (Reference Eqs. 2® 31* 2.32 

and 2.33) can be shown to approach -J- for Z * ^ . Also, 
for large arguments of the independent variable S and for 
“T * 4 i the following condition holds? 




Applying these two properties to Eqs. (2.29), (2. 30) 
and (2.35), it can be shown the inversion integrals for the 
radial, longitudinal and circumferential stresses fail to 
exist. 



CHAPTER III 


NUMERICAL EVALUATION OF INVERSE TRANSFORMS 

The numerical inversion of Eqs. (2,63) through (2.6b) 
is accomplished using a technique developed by Longman [ll] • 
Longman's method is based on Euler’s transformation of a 
convergent, alternating series [12]® A brief explanation of 
the method is given below. 

According to Euler's transformation, an infinite series 

\4 - V, + >4. - V s +■ V,f. - • • • ^ 

where 

\/^ >0 } J *44 sOj 

and 

- \C ^ -A * aW mh - a* } 

can be expressed as 

X 1C • *r vi - 4rA V. 4 4-dV, - * * - ( 3.1 ) 

Consider a function which oscillates about zero in 
such a way that the integral over each half-cycle is smaller 
in absolute value than that over the preceding half-cycle. 

The infinite integral of this function can be represented 
as an infinite alternating series where the ith term is 
the integral over the _ith half-cycle. Applying Euler's 
transformation to this series gives an accurate answer for 
a relative short interval of integration. The degree of 
accuracy is illustrated in Longman ' s paper [ll] , 
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To us® this method, it is necessary to determine the 
roots of the inversion integrands® The equations from which 
the roots are computed, based on asymtotic approximations 
of the Bessel functions, are derived in Appendix B, 

For large values of the radial parameter, successive 
magnitudes of the half-cycle integrations were found to 
initially increase. The largest magnitude existed for the 
third half-cycle. Using a series of seven terms ( integrals 
of the first seven half-cycles) to approximate the integral 
and applying Euler's transformation to the last five terms, 
gives for the value of the integral 

I — v; . v« ± £±v, *-rdV, + s* V) (3.2) 

In this form, each term retains its computed sign and the 
sign of the bracketed terms is the same as that of V $ . 

Plots of the inversion kernels , and , as func- 
tions of radial parameter are presented in Appendix C . 
Additionally, selected plots of various inversion integrands 
are presented. 

All numerical computations were accomplished by use 
of the Univac 1108. The method of integration was Simpson’s 
3/8 Rule. Routines for integration and evaluation of Bessel 
functions were obtained from the Univac library program, 

"MATH- PACK". 



CHAPTER IV 


RESULTS 

Figures 4 through 1? are plots of the nondimensional 
displacements and stresses as functions of the radial para- 
meter for three different values of the axial coordinate 2, 
Figures 4 through 7 are for the upper surface (Z * 0), 

Figs. 8 through 13 for the mid-plane ( Z * £** .5) and Figs. 

14 through 17 for the lower plane ( Z * A * jT*l • 0) • Each 
Figure has two curves j one corresponding to o< equal one and 
the other for o< equal two. The dimensionless parameter c< 
is defined as the ratio of hole radius to layer thickness. 

For all computations, the values V *.3 • 01 and^ *.020833 

were used. Figures 18 and 19 are sketches of the deformed 
state (heavy lines) superimposed on the undeformed state (light 
lines) for ©4 *1.0 and ©<*2.0, respectively. 

Figures 4, 8 and 14 are curves for nondimensional 
radial displacements, ^4 * Figures 4 and 8 are curves for 
tne upper surface and mid-plane, respectively. As shown, the 
radial displacements are smooth, monotonically decreasing 
functions. The curves in Fig, 14 give the radial displace- 
ment of the lower surface. The curves exhibit an initial 
increase , reach their maximum values at slightly greater 

than unity, and approach zero in the same smooth manner as 
the curves of Figs. 4 and 8. 
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The graphs for nondimensional longitudinal displace- 
ments, , are shown in Figs. 5 and 9 for the upper sur- 

face and mid-plane, respectively. Each curve shows negative 
values of displacement (upward) at T" • A . With increasing 
» each curve rapidly decreases in magnitude to zero, 
becomes positive and then approaches zero asymtotically. 

Figures 6 and 7 represent the nondimensional radial 
and circumferential stresses, and respectively, 

on the upper surface. It can be seen that the radial stresses 
are compressive and rapidly approach zero with 7 “/$, • The 
circumferential stresses are tensile and approach zero in 
an equally rapid manner. As one would expect, the radial 
stresses have the larger ' maximum values (at f* *A ) . 

Figures 10, 11, 12 and 13 are curves of the nondimen- 
sional radial, longitudinal, circumferential and shear stresses, 
respectively, on the mid-plane. Radial and circumferential 
stresses behave as they did on the upper surface and are 
approximately equal in magnitude. The longitudinal stress, 

CC^e- , is seen to rapidly decrease in magnitude. The stress 
is initially tensile, becomes compressive and then approaches 
zero. The curve for shear stress, g shows a rapid vari- 

ation in magnitude. The value of shear stress is initially 
zero , but quickly reaches a maximum and then approaches zero. 

Figures 15, 16 and 17 are curves for nondimensional 
radial, longitudinal and circumferential stresses, respectively. 
Each curve exhibits a singularity atiTsA. This feature will 
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be discussed later. The longitudinal stress is shown as 
initially tensile. A rapid change from tensile to compressive 
stress occurs with a subsequent tendency to aero. The curve 
of radial stress exhibits the same general properties. The 
curve for circumferential stress is monotonically decreasing 
and, with the exception of the singularity, possesses the 
same basic properties as the other curves of circumferential 
stresses. 

Figures 18 and 19 depict the deformed shapes of the 
layer for the two cases of « 1 and ©<.= 2, respectively. It 
can be seen that longitudinal displacement become insignificant 
with respect to radial displacement for increasing values of 

-r/tZ . 

The singularity that exists at T* * A, Z * j 4 can be best 
understood by considering the imposed boundary conditions. 

The shear stress on the cylindrical boundary is prescribed as 
zero. This fact, as shown by Eq. (2.7) » requires the following 
conditions 

On the lower surface, the longitudinal displacement, (jJ , is 
forced to be zero. This constraint on normal displacement 
gives rise to the following condition. 



From these two equations, it can be seen that at 
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X * inconsistent constraints result which will force a 
singularity. 

Curves for stresses and displacements are plotted 
using computed data over the range 1.05 < 3»5 • Below 

* 1.05» significant error was introduced because of trun- 
cation in the evaluation of the infinite integrals. The 
number of roots of the integrands that could be used in the 
integration technique described in Chapter III became limited 
due to their magnitudes. The roots became sufficiently large 
and caused computational overflow in the computer. As 
apnroaches £SL » the roots of each integrand increases without 
bound. At values of •f'/a close to unity, the first several 
roots take on very large values. Reference is made to Appen- 
dices B and C where the equations for computing the roots and 
curves showing integrand dependence on 7“ » respectively, 
are shown. With the exception of stresses at the singular point, 
all curves were extrapolated to r’/cL * 1.0« The extrapolation 
provided excellent results for radial displacements, 

Figs. 4, 8 and 14, where the data at f* A was known. Also, 
extrapolation of the curves for shear stress, , Fig. 13, 
provided good results at 7“ where the shear stress was im- 
posed as zero. 

It is interesting to note two features which are con- 
trary to what one might intuitively expect. From Figs. 15 and 
16, it can be seen that radial and longitudinal stresses are 
tensile at T" = A The longitudinal stress is tensile 
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because the surface of the cylindrical hole is negatively 
displaced (upward) but is constrained to be zero at the lower 
. surface. This extension gives rise to positive strain and, 
therefore, tensile stress. The cylindrical boundary being 
displaced upward is attributed to the absence of shear stress 
on its surface. The radial stress is tensile because the 
larger radial displacement imposed on the upper portion of 
the layer attempts by shear to displace the lower portion 
farther than the boundary condition will allow. This also 
explains the odd shape of the radial displacement curve, Fig. 
14, for the lower surface. 

Reasonable accuracy was obtained in checking the com- 
puted data by use of the following stress-strain relation 
for the case of axial symmetryi 

. 
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Figure 5® Nondimensional longitu- 
dinal displacement as a function 
of radial coordinate® 

f - 0 j8 * »Q1 4 > SB 9O2O833 
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Figure 6* Nondimensional radial 
stress as a function of radial 
coordinate* 


r .r o 


0 * .01 4 s *020833 



32 



Figure 7® Nondiraensional circum- 
ferential stress as a function 
of radial coordinate® 
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Figure 8* Nondimensional radial 
displacement as a function of 
radial coordinate* 


f - *5 fi * .01 * * « 0?0833 
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Figure 10, Nondimensional radial 
stress as a function of radial 
coordinate. 

f * .5 0 - .01 <f> = .020833 
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Figure 11* Nondimensional longi- 
tudinal stress as a function 
of radial coordinate* 

r « .5 0 « .01 $ * « 020833 
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Fieure 12. Nondimensional circum- 
ferential stress as a function 
of radial coordinate . 

f • ®5 0 * .01 * » ,020833 
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Figure 14® Nondimensional radial 
displacement as a function of 
radial coordinate® 

f *• 1 ® 0 0 » 
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Figure 16. Nondimensional longitu- 
dinal stress as a function of 
radial coordinate® 
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Figure 17. Nondimensional circum- 
ferential stress as a function 
of radial coordinate. 
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CHAPTER V 


SUMMARY AND CONCLUSIONS 

The problem of an elastic layer with a transverse 
cylindrical hole being subjected to an axially-symmetric 
radial deformation was investigated. The solution to this 
problem was achieved by use of extended Hankel transforms. 

The solution was based on transforming Navier’s equations of 
elasticity in the radial coordinate 7 " and solving the re- 
sulting ordinary differential equations in the axial coor- 
dinate z • The solution of these differential equations 
provided the displacement functions. The stresses were ob- 
tained by use of stress-strain relations. Inversion of the 
displacement and stress functions required the numerical eval- 
uation of infinite integrals. 

With the exception of the singularity, the curves for 
stresses and displacements were seen to behave in a very nor- 
mal manner. The imposed radial displacement on the cylin- 
drical boundary is obviously impractical. One could not ex- 
pect to maintain the sharp corner and remain within the confines 
of linear elasticity. The practical implication of this 
singularity (and the tensile radial stress) can be understood 
by considering a large disk of thickness 2^ with a transverse 
cylindrical hole of radius A « Consider this disk to be 
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shrink-fitted (frictionless) onto a shaft with a circum- 
ferential V-groove having the geometry shown in Fig* 3» It 
can be deduced from the results of this study that the zone 
of tensile radial stress would not be in contact with the 
shaft* Instead* the cylindrical surface of the disk would 
maintain a smooth shape in keeping with the boundary condi- 
tions of zero shear stress. At 1 = , the cylindrical sur- 

face would be perpendicular to the radial axis. 

It is further concluded that, while the forward trans- 
formations of Navier's equations and the analytical solutions 
thereof are relatively straightforward, the inversion pro- 
cess requires special attention. It is felt the technique 
presented herein is an acceptable method for numerically 
evaluating the infinite inversion integrals. 
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APPENDIX A 

PROPERTIES OF THE TRANSFORM 

The following properties of extended riankel trans- 
forms are taken from Scott and Miklowitz W* 


£ & fa r **-) ~ 



~1?r fat*) " 


4L 


- 



far, a ) Jr" * 



‘ £ *;a) 1~ STpJg, j - **>£ f*) 
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(A3) 

(A4) 
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Additional properties ( see, for example. Reference 

2) ares 



2 

(A‘;) 


O 

(A6) 

= 

-5-^/ 

(A?) 



(A8) 

>*- [£>, = 
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Par las [13] showed by contour integration that the 

I 

following condition is true* 
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APPENDIX B 

ROOTS OF THE INTEGRANDS IN THE INVERSE TRANSFORMS 

In order to find the roots of the inverse integrands , 



S Js 

(Bl) 


JiV* **) •*- 



(B2) 

it is only necessary to determine the roots of and 


which are 

defined as followss 


/2 f *jm) 

- j# Ar) y f&®) - J, {$&) y$ ($r>) 

(B3) 

fa*;*) 

- J, tsr) y, fa 1 ) - J, fat) y, far) 

(B4) 


The derivation of the expressions for the roots 
utilizes Stoke ’ s method (see, for example. Reference 14)® 

The following are asymtotic approximations for Bessel func- 
tions of the first and the second kind [l 5 ] 1 

Pft*) c&s-y* -Qfxid (B5) 

#** *»®n> 

Ytor> ^ y^ [pftx) j ^ 


(B6) 
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where 

Qftx) ^ 


z ........ 

z/ ^f/) 8 y/ ^0r)f 

y?*7 _ &?*/X?r-*)L 

m ' ~ 3/ y*x ) 5 






Roots of Rofatrall. 

Given Eq. (B3) and letting ’fio « svz. and y® rfa. } 
Eq. (B3) becomes? 

^5 y) - J. /««) y, fa) - -V A> y> O'*) . 

It follows that Eqs. ( B5) » ( B6 ) , ( B7 ) » ( B8) and (B9) can be 
written? 

Jfi fo) “ *’ $(*> t 



(By) 
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(BQ) 

(BIO) 
(Bll) 
( B12 ) 



54 







where 




where 



Setting £q» (BIO) equal to zero and realizing Lhe 
equalities that exist with the s/gi* % “** J ^ 
terms , the following simplification results » 


( f. -Iffi - % Tf e 




For positive values of -ft, » the zeros of Eq, (B15) 
are represented bys 



J&js ^ i *, %j, • - - . 


Rearranging to solve for $ in terms of ^ gives? 
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where % represents the nth root of Eq. (BIO), 
Letting 




£ + 


£ 


£& 


and assuming 

^ ‘ " r i* 

it can be shown the expression for computing the nth root 
of is * 


(Ble; 
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Hoots of 

Using the same procedure as before, the expression 
for computing the nth root of £(*,**) is* 


rffa*) 


r-/ $Y Saht} 

V*yj^g/ 

/*<# Jr* * JP*sU3 


Comment . 

The expression for the roots of J2,(syfjd.^ , Eq. (B21 ) , 
is verified by Gray and Mathews [16] « Data generated from 
Eq® ( B20) have been verified by the curves of the functions 
presented in Appendix C. There is slight error in the first 
root due to the asymtotic approximations of the the Bessel 
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functions being; for large arguments. This error* however* 
had no effect on the results presented herein. 
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APPENDIX C 


CURVES FOR INVERSION KERNELS AND SELECTED INTEGRANDS 


Curves for the inversion kernels K and K, and 
inversion integrands tXK s and ^1 K® are presented as func- 
tions of the independent variable S. These curves are 
representative in nature and are not intended to imply 
specific importance. 

The terms W and <3jt are defined in Chapter II of 
the text and K 0 and are defined asi 


K„ 

Ki 

where 
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Figure Cl « Zero-Order inversion 

kernel versus independent variable , 
r/a- 1 » 05 * 


Figure C2. Zero-Order inversion 

kernel versus independent variable* 
r/a^l , 25® 








Figure C3« Zero-Order inversion 

kernel versus independent variable, 
r / a™ 1 ® 5 ® 



61 



Fieoire G4« Zero~Order inversion 

kernel versus independent variable* 
r/as2»5» 



s 


Figure C5® Zero-Order inversion 

kernel versus independent variable * 

r/a=3« 5® 
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Figure C6a First-Order inversion 

kernel versus independent variable* 
r/a=l « 05 
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Figure C7, First-Order inversion 

kernel versus independent variable* 
r / a— 1 ® 2 5 * 



Figure C8® First-Order inversion 

kernel versus independent variable* 
r/as4®5» 
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Figure C9® First-Order inversion 

kernel versus independent variable » 
r/a=2» 5® 
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Figure Cll» Inversion integrand versus 
independent variable* radial dis- 
placement, r/a=l * 1 * 





Figure Cl 2* Inversion integrand versus 
independent variable, radial dis- 
placement, r/a*3.5„ 
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Figure Cl 3® Inversion integrand versus 
independent variable, longitudinal 
displacement, r/a=l ,1, 
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§ 


Figure C14. Inversion integrand versus 
independent variable, longitudinal 
displacement, r/a=3* 5« 
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